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Abstract 

Rapid rebinding of molecular interaction partners that are in close 
proximity after dissociation leads to a dissociation and association ki¬ 
netics that can profoundly differ from predictions based on bulk re¬ 
action models. The cause of this effect can be traced back to the 
non-Markovian character of the ligand’s rebinding time probability 
density function, reflecting the fact that, for a certain time span, the 
ligand still ’remembers’ the receptor it was bound to previously. In this 
manuscript, we explore the consequences of the hypothesis that initial 
binding and consecutive rebinding give rise to a bond lifetime density 
that is non-Markovian as well. We study the combined effect of the two 
non-Markovian waiting time probability densities and show that even 
for very short times the decay of the fraction of occupied receptors devi¬ 
ates from an exponential. For long times, dissociation is slower than an 
exponential and the fate of the the steady-state bound receptor fraction 
critically depends on the extent of the deviation, relative to the rebind¬ 
ing time density, from the Markovian limit: The population of occupied 
receptors may either decay completely or assume a non-vanishing value 
for small and strong deviations, respectively. Furthermore, we point 
out the important role played by fractional calculus and demonstrate 
that the short- and long-time dynamics of the occupied receptors can 
be naturally expressed as well as easily obtained in terms of fractional 
differential equations involving the Riemann-Liouville derivative. Our 
analysis shows that cells may exploit receptor memory as mechanism to 
dynamically widen the range of potential response-patterns to a given 
signal. 
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1 Introduction 


Growing evidence [H HI [25l [18] suggests that a ligand’s diffusive motion 
results in a receptor-ligand interaction pattern and a network response 
further downstream that can deviate substantially from the picture 
based on deterministic first-order kinetics. A major cause of these 
differences has been identified as the non-Markovian rebinding time 
probability density function (pdf) that arises in stochastic, spatially- 
resolved accounts of the ligand’s motion [25l|T8]. This rebinding time 
pdf may be considered as a ’deformed’ version of the memoryless ex¬ 
ponential waiting time pdf that corresponds to the first order kinetics 
of a deterministic mean-field description. In fact, the rebinding time 
pdf iprebit) for 1-dimensional diffusive motion in the presence of a par¬ 
tially absorbing boundary takes in the Laplace domain the form (see 
appendix A) 


/■oo 1/2 

>C[V'reb(t)](s) := Aehis) := / V’reb(t)e“'’*dt = , ^ --, (1) 

Jo Kj^ + S^/^ 

where ku = ^ parameter whose strength is determined by the 

ratio of the intrinsic association constant Ha to the diffusion constant 
D. Note that kd has the dimension of an inverse time. 

By contrast, the exponential waiting time pdf that is commonly 
employed to describe the dissociation reaction (and association as well 
in a non-spatial deterministic approach), 

iliit) = ( 2 ) 

reads in the Laplace domain as 


where Hd refers to the intrinsic dissociation constant. Both density 
functions (Eqs. ([T]), dSj)) are representatives of a family of pdf that are 
parametrized by a parameter 0 < cr < 1 muoiEi] 


^(s) 


-|- S° 


(4) 


Any PDF with cr < 1 describes memory effects: Because a ligand, that 
just dissociated from a receptor, will remain close to that receptor for a 
certain time span it is more likely to bind the same receptor again. The 
spatial correlations between receptor and dissociated ligand translate 
to a time evolution of the ligand’s (survival) probability that is non¬ 
local in time and hence give rise to a dynamics that can be naturally 
cast as a fractional differential equation. Fractional calculus [T9| [17] 


2 





plays a prominent role in anomalous diffusion and non-equilibrium sta¬ 
tistical mechanics of complex systems but it also 

naturally arises in the context of normal diffusion. To see this in more 
detail, let us recall that the rebinding PDF '0reb(i) of a single ligand 
is related to the survival probability of a molecule that undergoes a 
1-dimensional diffusive motion in the presence of a partially absorbing 
boundary condition (BC) at zq = 0 by (cf. Appendix A) 


'0reb(^) — 


9S'rad(i|-20 = 0) 
dt 


(5) 


Henceforth, we will use the notation S{t) := S'rad(i|-2o = 0) for conve¬ 
nience. From Eqs. o, ® and (HSl) we can conclude that 


1 — sS{s) = 




This relation reads in the time domain as 

_ _ l/2^1-l/2 


dt 


Sit), 


( 6 ) 


(7) 


where DI ° denotes the Riemann-Liouville (R-L) fractional differential 
operator that is defined as [mmni] 

1 d I 


:= := 


/M 


r(o-)5tyo (t-T) 


1—(T 


dr, 0 < cr < 1. 


( 8 ) 

The dehnition of the R-L derivative (Eq. ®) makes it evident that 
the time evolution of S{t) (Eq. (O) is indeed governed by an equation 
that is non-local in time. Eq. o generalizes the well-known Markovian 
time evolution 

^ = -^DS{t) (9) 

that leads to a Poissonian form of the rebinding pdf. Solutions of 
evolution equations that take the form of a fractional relaxation law 

[Timiin] 


dt 




( 10 ) 


are given in terms of Mittag-LefHer (M-L) functions O [13] $(f) = 
Ea[—{Kt)‘^]- M-L functions may be thought of a generalized exponen¬ 
tials that are defined by the series expansion 


Ec{-x) = 


{-xY 


n—0 


T{na + 1) ’ 


( 11 ) 


which makes it evident that for a = 1 one recovers an exponential 
Ei{—t) = exp(—t). The Laplace transform is 


C 


Eai-iKt)'^) 


1 


s s'^ + k'^ 


( 12 ) 
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and hence it follows that the pdf Eq. 0 is the Laplace transform 
of the negative time derivative of the Mittag-LefHer (M-L) function 

Erj[-{KtY] [TO]: 


V'(i) 


dEY-[K,tY] 

dt 


(13) 


Because Ei/ 2 [—= exp(t)erfc(-\/t) we immediately obtain from 
Eq. ([7|) the known expression for S{t), cf. Appendix A, Eq. (ITOll . 
For small arguments, the M-L function behaves like a Kohlrausch- 
Williams-Watts (KWW) function, i.e. a stretched exponential ~ 

exp(— [M] [TO] (cf. Eq. (fTTIl l. while for large arguments it decays 
like an inverse power law. More precisely, one has [TO] HI] 


Ej-n = 


t 


r(l-a) r(l-2a) 


t oo. (14) 


We have seen that fractional evolution equations (Eq. ([7]) ) that re¬ 
flect deviations from a Markovian waiting time pdf naturally appear 
in the decription of ligand rebinding, even in the case of normal diffu¬ 
sion. However, the reverse process, dissociation is usually taken into 
account by a memoryless pdf for the bond lifetime, even in particle- 
based stochastic approaches (however, note, for instance, Refs. [71 ITT]'). 
Here, we suggest to abandon the assumption of the Markovian prop¬ 
erty. A deviation from a memoryless bond lifetime pdf is motivated 
by several lines of experimental evidence suggesting that the receptor 
state does not remain unaltered upon initial rebinding, for instance 
due to receptor deformation or conformational changes, the formation 
of signaling intermediates, involvement of the actin skeleton and clus¬ 
ter formation and other mechanisms [UHIIIIT]. In the picture that 
we will consider, the bound state undergoes a diffusion-like ’motion’ in 
an abstract internal ’space’ whose properties are dynamically defined 
by the local biochemistry near the receptor. We will show that the 
interplay of these two different diffusive processes alters substantially 
the apparent dissociation kinetics and can even lead to behavior that 
has no counterpart in models that only take into account either the 
non-Markovian character of the ligand rebinding or the bond lifetime 
alone. 

Let us recall how an altered receptor state may lead to a bond 
lifetime pdf that is sufficiently well described by Eq. (J?]) [^ [7] . To this 
end, one assumes that the receptor modifications imply a bond lifetime 
pdf that is given by a superposition of exponentials 


One further assume that the sum can be approximated by an integral 



( 16 ) 
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Then, using mild assumptions, one can derive (cf. Appendix B) that, 
asymptotically, the pdf behaves like a power-law 

Employing Tauberian theorems liHU one can show that 

V;(s) ~ 1 - (rs)'", (18) 

where r denotes a constant having dimension of time. Now, for large 
times, i.e. small s or for small t we finally may rewrite Eq. (HI as 
Eq. O, identifying r ^ = k. 

Therefore, we will henceforth assume that the bond lifetime pdf is 
indeed given by Eq. We consider the following system: Ligands 
diffuse in a 3 dimensional space that is semi-infinite, bounded by a two- 
dimensional plane homogeneously covered with receptors. We adopt 
the self-consistent stochastic mean-field theory approach of Ref. [5] that 
allows to track the motion of the ligand accurately, while the receptors 
are described in a non-spatial fashion by the fraction of occupied recep¬ 
tors. We point out that the formalism is flexible enough to include also 
a spatially-resolved description of the receptors [9]. Here, however, we 
aim at emphasizing the consequences of a non-Markovian bond lifetime 
density and its interplay with the rebinding pdf and therefore limit our 
analyses to a non-spatial receptor description. 


2 Theory 

First, we focus on the Markovian case [8]. We consider a membrane 
that extends infinitely in the x — y plane {z = 0) and that is uniformly 
covered with receptors. We assume that all ligands are initially bound. 
The fraction of occupied receptors is denoted by p{t) = TZ{t)/TZo where 
Ti{t) and TZq refer to the density of bound receptors at time t and the 
constant mean receptor surface density per unit area, respectively. Its 
time evolution is governed by [8] 

= -kospit) + konp{0,t)[l -p{t)], (19) 

where p{z,t) describes the concentration of ligands at a distant z from 
the membrane. To solve Eq. (HU) one has to express the ligand density 
p(0,t) and hence the rebinding rate 

7(i) = konP(0, t) [1 - p{t)\ (20) 

in terms of p{t). To this end, we follow the self-consistent mean-field 
theory approach described in Refs. [8] and assume that p(0) <C 1 and 
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hence p{t) <C 1 for all times t so that the rebinding rate reduces to 
7 (t) = fconp(0,t). We now aim at expressing ^{t) in an alternative 
form that incorporates spatially-resolved information about the diffu¬ 
sive motion of the ligand. The fraction of occupied receptors which is 
converted to the unbound state between time r and r -|- dr is given by 
koSp{T)dT. A dissociated ligand begins to undergo an effectively one¬ 
dimensional diffusive motion that starts at the mebrane z = 0 between 
time r and r -I- dr. Let — r|z = 0) denote the probability that a 
ligand is still unbound when the time span t — t has elapsed since the 
dissociation. Then, the rate of absorption is given by by the rebinding 
time pdf (see Appendix C) 

'iprehit) := R^e,dit - t\zo = 0) - — ( 21 ) 

and it follows that the rebinding rate can be written as convolution 

lit) = ^off [ dTp{T)Rj.g^^{t - r|0). (22) 

Note that, despite the formal resemblance of Eqs. m and m , we 
emphasize that the context of both relations is quite different, because 
here we consider a many-molecule system, instead of a single ligand- 
receptor pair. Also, the self-consistent mean-field theory approach is 
more flexible and can be applied to more general cases [9]. 

As it stands, Eq. dm) describes memoryless dissociation based 
on the exponential waiting time pdf Eq. ©■ To introduce a non- 
Markovian intrinsic bond lifetime, we perform the substitution 

kos ^ (23) 

1 - V'(s) 

in the Laplace domain, where ip{s) is given by Eq. (|3]) with k = /coff- 
Clearly, for ct = 1, one recovers s'!/'(s)/[l — V'('S)] = feoff- The substitu¬ 
tion rule Eq. (1251) has been utilized before mm- However, Refs. mm 
have not identified the action of the operation of Eq. (l23l) in the time 
domain. Using [l4] 

= (24) 

it turns out that it corresponds to 

- feoffP(f) — -k^gDl~'^p{t) (25) 

where Dl~^p{t) refers to the R-L derivative (Eq. ([5])). Now, combining 

Eqs. (fT2]) . (1121) . (1221) . (1251) . (1251) and (1521) we obtain 

p{s) = s-^piO) - -ITIpW, (26) 
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or, equivalently, 


p{t) = p{0) - k^ffDl 


El/2[-{KDt)^^'^] *p{t) 


(27) 


where * denotes a convolution integral. 

Eq. ^ represents a natural generalization of a rate equation ap¬ 
proach, abandoning the assumptions of both a memoryless rebinding 
and dissociation waiting time pdf. In fact, considering the limit ct —^ 1 
and substituting Ei/ 2 [—{KDty^^] = exp(—we 

obtain from Eq. (1271) the deterministic Markovian mean-field result 


p{t) = p{0) 


^-(K.n+k„„)t ^off 

KD + fcoff 


Kp 

KP + fcoff 


(28) 


Next, we turn to the general case and focus on the short and long time 
limit. We will show that in most cases the dynamics is described by 
the same fractional relaxation law (Eq. (jlOp l. only the order of the R-L 
derivative changes from case to case. Hence, as a benefit, the solutions 
can be virtualy read off, a fact which underlines the power and elegance 
of the chosen formalism. 


3 Results 


3.1 Solution for short times t —)■ 0 

For very short times, the effect of the rebinding should be negligible. 
This expectation is well reflected by Eq. 1271); Because short times 
imply that Ei/ 2 [—{Kpty/^] —>• 1, it follows that the equation of motion 
assumes the form of a fractional relaxation law, 

^ = -k^sDl-Xt), (29) 

for which we can swiftly read off the solution (cf. Eqs. (flOl) and (fTTI) ') 


P(.t) 


p{0)E„[-{ko«ty 

p(0)fl- 


r(i + a) 


jkosty'^ 
r{l + 2a) 


(30) 

fort —0. (31) 


Note that to arrive at Eq. IMD, we invoked the commutation relation 
,d/dt]f{t) = 0 that does hold if /(O) = 0. Since in the case 
considered here we have /(<) = J*p{T)dT, the condition is fullfilled. 
As we know from the behavior of the M-L function (and can see from 
Eq. ini)), for short times the fraction of occupied receptors decays like 
a stretched exponential pit) ^ p(0) exp[—(fcoffi)'^]- 
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We emphasize that the solution given by Eq. o is not only 
valid for short times, but in addition in parameter regimes in which 
= TZokon/^/D becomes small, i.e. when either the receptor density 
is small or when rebinding effects are weak, i.e. the ratio of on-rate to 
diffusion constant becomes small. Because for ct = 1 the correspond¬ 
ing M-L function in Eq. (l30l) is a decaying exponential £’i[—(fcofft)^] = 
exp(—fcofft), we recover the result from [5] for short times. However, 
here we see that abandoning the assumption of a memoryless bond 
lifetime leads to an decay that deviates from an exponential even for 
short times or for small kd- 


3.2 Solutions for long times t ^ oo 


Now we will show that the combined effect of rebinding and receptor 
memory gives rise to behavior that has no counterpart in cases in which 
one uses only one of the two non-Markovian densities. The long time 
limit, or equivalently s —>• 0, implies that (cf. Eq. (fl^ l 

c(E,/,[-iKDty/^]] (32) 

As a result, Eq. (1271) again takes the form of a fractional relaxation 
law, as in the case of short times, but the order of the R-L derivative 
has changed (cf. with Eq. ((2^ 1 


dp{t) 

dt 


A/2 

''D 


j^l-a+1/2 


p{t). 


(33) 


Here, we again emphasize that the found relaxation law is not only 
valid for long times, but also for regimes in which kd becomes large, 
i.e. where rebinding effects becomes strong. Furthermore, the form 
of Eq. (1551) makes it evident that there is now a crucial difference 
whether cr < 1/2 or cr > 1/2 i.e. whether the bond lifetime is to a 
larger extent non-Markovian than the rebinding time pdf or not. In this 
sense, diffusion sets the scale for what should be considered a strong or 
weak deviation from memoryless behavior. Note that this scale is not 
absolute: If the ligand undergoes anomalous diffusion HHHaESlEl], 
the scale will not be 1/2 anymore. 

We now consider several special cases separately. First, let us focus 
on the Markovian limit cr = I. Then, we immediately obtain from 
Eq. dMI) 


P{t) = Ei/2 




(34) 


which is the result given in Ref. 019], as expected. In particular, we 
observe a power law decay for long times p{t —> oo) ^ \j\piri. 
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Next, we analyze the case 1 > cr > 1/2. It follows that ct — 1/2 > 0 
and hence we can read off again immediately the solution from Eq. (1331) 


P{t) 


P{Q)Ecr-l/2 

p{0) 


, 1/2 




-{a- l/ 2 )]t-i /2 


Kd 

kl-T[l - 2{a - l/2)]t2-i 


where we employed the expansion of the M-L function for large argu¬ 
ments Eq. da. We note that the fraction of occupied receptors will 
completely decay for long times, as in the Markovian case cr = 1. How¬ 
ever, it can also be seen from Eq. (1551) that for general 1/2 < cr < 1 
the decay can become arbitrarily slow for cr —>■ 1/2 and hence cr = 1/2 
seems to play a special role that we are going to consider now. 

As previously, the asymptotic solution p(t —>■ oo) for the limiting 
case cr = 1/2 can easily be obtained from the general relaxation law 
Eq. (|^ and using Eq[—x) = 1/(1 -I- a;) 


( 35 ) 


p(f-» oo) =p(0) , 1 / 2 - (36) 

Thus, we find that the steady-state receptor occupancy is not vanishing 
in the long time limit, in stark contrast to the case of an exponential 
decay of the intrinsic bond lifetime. We point out that Agmon and 
Weiss found a corresponding behavior for the case of a single pair in 
one dimension [2]. But as mentioned before, the formalism adopted 
here is more general Moreover, we establish a connection to 

3D/2D membrane binding as well as fractional calculus and the crucial 
role of the fractional relaxation law. 

Finally, we turn to the case cr < 1/2. This is the only case for 
which we cannot read off the solution immediately. The reason is that 
the nature of the R-L derivative changes qualitatively. In fact, let be 
a = 1/2 — a, then a > 0 and one has [H] 


d4 


pit) = 


52 


p(r) 


r(l — a) df^ Jq [t — t)' 


-dr. 


(37) 


Therefore, we Laplace transform Eq. (1331) (or, alternatively, start from 
(l26l) l and obtain 


Pis) 


A/2 


p(0) 


„< 7 — 1 




'-off' 


A/2 


p(0) 


1 -f K 


1/2 

D 


1 


(38) 
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The corresponding expression in the time domain follows from Eq. (HH) 


pit) =p(0)|l - Ei/2-a 
In particular, we hnd that 




lim p{t) = p(0) 

t—¥00 

the initial fraction of occupied receptors is restored completely. 


(39) 


(40) 


Appendix A 

We consider a molecule undergoing 1-dimensional diffusive motion in 
the presence of a partially absorbing BC. The pdf grad{z,t\zo) that 
yields the probability gra.d{z,t\zo)dz to find the molecule between z 
and z + dz at time t, provided that it was initially located at zq at time 
to = 0, is the GF of the diffusion equation 


dgrs,d{z,t\zo) _ j^d‘^gTad{z,t\zo) 


dt 

subject to the radiation BC [31 Eg 
,%rad(^,t|2o) 




(41) 


D- 


dz 


=0 = Ka3rad(z = 0,t|zo)- 


(42) 


The GF takes in the Laplace domain the form |3l p. 358] 

q-\/s/D\z-zo\ ^-^JTpD{z+zo) ^ g-\/ s/D(z-za) 

5rad(^,s|zo) = 


2^/ Ds 2'J Ds 

The survival probability 5'rad(t|-Zo) is defined as [! 

POO 

5'rad(iko) = / grad{Z:t\zo)d. 

Jo 


D^/s[^/s + K,a/ " n /^] 
(43) 


and hence one hnds 


*^rad(s|zo — 0) 


,1/2 


(44) 


(45) 


S sl/2 + k^/2 ’ 

where kd = Eq. (1451) becomes in the time domain [U Eq. (29.3.43)] 

-S'rad(t|2o = 0) = exp(K_Dt)erfc(VKDt). (46) 

Furthermore, from Eqs. dUD and (|42ll one may derive the relation 

^^rad (t 12o ) 


dt 


= Hagri,d{z = 0,tlzo)- 


(47) 
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Finally, using Eqs. 63]) and dUl) we obtain 


V^reb('5) — 


55rad(i|zo = 0) 


dt 


. 1/2 


+ sl/2 


(48) 


Appendix B 

In this Appendix we recall Enzma how the form of the density 
that appears in Eq. can be derived from rather mild assumptions. 
Eirst, one assumes that a typical rate k is characterized by an Arrhenius 
factor 

K oc e“'^, (49) 

where fc, T refer to the Boltzmann constant and the absolute temper¬ 
ature, respectively. The activation energies E themselves are assumed 
to be distributed according to an Arrhenius factor too 

(j){E) oc e~^, (50) 


where Eq denotes a energy reference scale. Using Eq. (P|. one may 
express B hy and substitute the obtained expression in Eq. (1501) to 
hnd 

(/>(k) oc (51) 


where a := kT/E q. Eq. (1161) defines a Laplace transform where k plays 
the role of the time variable. Hence, knowing Eq. m, one may employ 
£[t“-i] = r(a)/s“, a > 0 [TJ Eq. (29.3.7)] to finally arrive at Eq. (ITTl) . 


Appendix C 


We consider first a single ligand-receptor pair. The rate of recombina¬ 
tion i?rad(/|2:o) = —9Sra.d{t\zQ) / dt Can be written as convolution [00] 


^rad(/bo) = 0-Rabs(/|2:o) + (l-0) [ ^abs(T|Zo)i?rad(i-T| A)dT for A 0 

(52) 

where i?abs(/|2:o) refers to the rate of first encounters (first-passage time 
pdf) 

R.Ut\zo)=-Bt^, (53) 

at 


and S'abs(/|2o) denotes the survival probability in the presence of an 
absorbing BC. Furthermore, A is a microscopic length scale and 0 
denotes the probability that a single molecule pair reaction takes place 
when they encounter. However, here we are interested in a situation 
where we consider many receptors. Eq. dsa) remains valid O |9] , but 
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the probability 0 has to be adapted properly, 
we define 

TZokon^ 


Following Refs. 019], 
(54) 


Next, we Laplace transform Eq. (|5^ to arrive at 


.Rrad(sl^o) = 0.Rabs(sko) + (1 “ 0).Rabs (s |^ 0 )^rad (s| A) . (55) 

Setting zq = A, we obtain from Eq. (1551) an expression for i?rad(s|'2^o = 
A) and substite it back in Eq. (1551) to find 


.^rad('5|'2^o) — 


0.^abs(5|-2'o) 


1 - (1 - 0)Rabs(s|A) 


(56) 


Using Ra_hsis\zo) = e [5] and Taylor-expanding up to first order 

in A yields 

QRs,hs{s\zo) 


^”d(®l^°)"0 + (l-0)A/v^sl/2 
Einally, we take the limit A —^ 0 (cf. Eq. (l54ll l 

hm ^ ^ 


A—fO 0 ^on"^0 


(57) 


(58) 


to arrive at 


V^reb('5) — Rrad('5|^0 — 9) — —1/2 / ’ '— ”^0 ^on j D. 

1 + Kjj 

. 

In particular, this result means that the rate of recombination is the 
negative time derivative of the M-L function Ei/ 2 [—TT-ofeon/V^t^^^]- 
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